Non-Markovian quantum dynamics: Correlated projection superoperators and 

Hilbert space averaging 
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The time-convolutionless (TCL) projection operator technique allows a systematic analysis of the 
non-Markovian quantum dynamics of open systems. We introduce a class of projection superopera- 
tors which project the states of the total system onto certain correlated system-environment states. 
It is shown that the application of the TCL technique to this class of correlated superoperators 
enables the non-perturbative treatment of the dynamics of system-environment models for which 
the standard approach fails in any finite order of the coupling strength. We demonstrate further 
that the correlated superoperators correspond to the idea of a best guess of conditional quantum 
expectations which is determined by a suitable Hilbert space average. The general approach is 
illustrated by means of the model of a spin which interacts through randomly distributed couplings 
with a finite reservoir consisting of two energy bands. Extensive numerical simulations of the full 
Schrodinger equation of the model reveal the power and efficiency of the method. 
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I. INTRODUCTION 

Realistic quantum mechanical systems are influenced 
through the coupling to an environment which contains 
a large number of mostly uncontrollable degrees of free- 
dom. The unavoidable interaction of an open quantum 
systems with its environment gives rise to the mecha- 
nisms of damping and dissipation, and to a strong and 
often rapid loss of quantum coherence. Applications of 
the theory of open quantum systems p] are found in al- 
most all areas of physics, ranging from quantum optics 
to condensed matter physics Jj and chemical physics 
from quantum information Hjto spintronics §. More- 
over, the theory of open quantum systems provides the 
foundations of quantum measurement theory Q, deco- 
herence @ and the emergence of thermodynamic behav- 
ior 0. 

In a microscopic approach one regards the total sys- 
tem, which is composed of the open system S and its 
environment E, as a closed quantum system following 
a Hamiltonian time evolution. One of the central goals 
of the theoretical treatment is then the analysis of the 
dynamical behavior of the populations and coherences 
which are given by the elements of the reduced density 
matrix ps(t) — tr Epif)- Here, p(t) denotes the density 
matrix of the composite system and tr e the partial trace 
taken over the environment. 

In the Markovian regime a complete mathematical the- 
ory is available which is based on the concepts of com- 
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pletely positive quantum dynamical semigroups and the 
corresponding Markovian master equations in Lindblad 
form 0EII1I3- However, if the physical conditions 
underlying the Markov approximation are violated one 
has to cope with strong non-perturbative and memory 
effects and the theoretical and mathematical treatment 
of the reduced system dynamics is typically much more 
involved. 

A systematic approach to non-Markovian dynamics is 
provided by the projection operator techniques [TEl ITt^ 
which are extensively used in nonequilibrium thermody- 
namics and statistical mechanics Il7l 0. The key con- 
cept of these techniques consists in the introduction of 
a certain projection superoperator V which acts on the 
operators of the state space of the total system. The su- 
peroperator V formalizes the idea of the elimination of 
degrees of freedom from the complete description of the 
states of the total system. Thus, if p is the density ma- 
trix of the composite system, the projection Vp serves 
to represent a simplified effective description through a 
reduced set of variables. For this reason, the projection 
Vp is called the relevant part of the total density ma- 
trix, while the complementary projection Qp — p — Vp is 
referred to as irrelevant part. 

With the help of the projection operator techniques 
one derives closed equations of motion for the relevant 
part Vp(t) from which one obtains approximate master 
equations by means of a systematic perturbation expan- 
sion with respect to the system-environment coupling. 
We shall concentrate in this paper on a specific variant of 
the technique which is known as the time-convolutionless 
(TCL) projection operator method [H HJ il |H |23| . 
The advantage of this formulation is that it leads to dy- 
namic equations for V pit) which are local in time and in- 
volve an explicitly time-dependent generator. A general 
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account of projection operator methods and, in partic- 
ular, of the TCL approach and its applications may be 
found in Ref. 

In the standard approach to the dynamics of open sys- 
tems one chooses a projection superoperator which is de- 
fined by the expression Vp = ps® Pe, where pe is some 
fixed environmental state. This superoperator projects 
the total state p onto an un-correlated tensor product 
state. Since pe is considered as fixed, it implies that the 
elements of the reduced density matrix ps(t) represent 
the relevant variables used for an effective description of 
the reduced system dynamics. This ansatz for the pro- 
jection superoperator V is widely used in studies of open 
quantum systems. It has been used to derive Markovian 
and non-Markovian quantum maste r eq uations for many 
applications (see, e. g., |5E |H il Moreover, 
non-Markovian generalized master equations for the re- 
duced density matrix have been developed on the basis 
of phenomenological considerations 29]. 

The paradigm of these and many other approaches is 
the usage of the reduced density matrix ps(t) as dynami- 
cal variable for which appropriate (exact or approximate) 
dynamic equations are to be developed. However, it is 
important to realize that the projection operator tech- 
niques are much more general and flexible, and that they 
offer many further possibilities for the construction of 
suitable projection superoperators. The only formal con- 
dition which must be satisfied in order to apply the tech- 
niques is that V is a map which operates on the total 
state space and has the property of a projection opera- 
tor, i. e., V 2 = V. This is a very general condition which 
can be fulfilled in many different ways. 

We mention two examples: In the analysis of classi- 
cal stochastic processes one considers a projection of the 
form VX(t) = (X(t)) which takes any stochastic pro- 
cess X(t) to its average (X(t)) With this choice 
the TCL technique leads to a cumulant expansion of the 
dynamic equation for the average [8fll IsiL l82l Isflj . In 
nonequilibrium thermodynamics a further projection of 
the form Vp = Pdiag is often introduced which maps any 
density matrix p to its diagonal part paiag i n & suitably 
chosen basis, yielding the famous Pauli master equation 
in lowest order of the TCL expansion . 

In the present paper we shall construct a class of 
projection superoperators V which enable the non- 
perturbative treatment of highly non-Markovian pro- 
cesses in open quantum systems. These superoperators 
project the state of the total system onto a correlated 
system-environment state, i. e., onto a state which con- 
tains statistical correlations between certain system and 
environment states. Thus, we give up the paradigm of us- 
ing the reduced density matrix as the dynamical variable, 
and enlarge the set of relevant variables to account for 
statistical correlations which are responsible for strong 
non-Markovian effects. The idea of introducing addi- 
tional variables has been realized in different ways and 
used in various contexts [2 OS EE EE ISE Here, we 
implement this idea directly in the definition of a projec- 



tion superoperator and connect it with the method of the 
TCL technique. This connection enables us to determine 
higher order corrections in a systematic way and, hence, 
to assess the quality of the approximations obtained. 

Recently, an entirely different approach has been sug- 
gested, the Hilbert space average method (HAM) |9]lirJ. 
This method employs the concept of a best guess for 
conditional quantum expectation values. It is based on 
the determination of a conditional Hilbert space aver- 
age. The method provides us with a systematic principle 
to estimate quantum expectation values conditioned on 
prescribed values for the expectations of a given set of 
operators on the total state space. HAM can be used to 
construct effective equations of motion for the given set 
of operators and, hence, yields an alternative approach 
to non-Markovian quantum dynamics. It will be shown 
here that the method of the Hilbert space average and the 
projection operator techniques which are based on the 
class of correlated projection superoperators are closely 
related. In fact, we are going to demonstrate that HAM 
represents the lowest order of the TCL expansion corre- 
sponding to this class of superoperators. 

The application and the efficiency of our approach will 
be illustrated and discussed here by means of a specific 
system-reservoir model. The model consists of a spin 
which interacts with two environmental energy bands 
through a set of random couplings 0, Ell • This model 
exhibits an unexpected feature. Namely, it turns out 
that the usual Born-Markov approximation fails for this 
model, although the standard Markov condition is satis- 
fied, i. e., although the width of the environmental two- 
point correlation function is small compared to the re- 
laxation time. By contrast, it will be demonstrated by 
means of a comparison with the numerical solution of 
the full Schrodinger equation of the model that the TCL 
expansion using the correlated projection superoperator 
yields accurate results already in lowest order of the per- 
turbation expansion. 

The paper is organized as follows. The class of cor- 
related projections will be introduced in Sec. [H] This 
section also contains a brief general account of the pro- 
jection operator techniques, and introduces the basic dy- 
namic equations as well as the perturbation expansion 
of the master equations for the relevant variables. The 
principles and equations of the Hilbert space averaging 
method are outlined in Sec. II I II where we will also de- 
scribe the connection between the principles of HAM and 
the structure of the correlated projection superoperators. 
Sec. IIVI contains the application of the general concepts 
developed here to a specific system-reservoir model. We 
shall discuss in detail the origin of the failure of the Born- 
Markov approximation, and show that and why the new 
projection superoperators yield an efficient and accurate 
approximation of the dynamics. Finally, we draw our 
conclusions in Sec. 
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II. PROJECTION OPERATOR TECHNIQUES 
A. Projection superoperators 

We consider an open quantum system S with state 
space Tis which is coupled to an environment E with 
state space He- The Hilbert space of states of the 
composite system is given by the tensor product H = 
Hs ®He- We assume that the dynamics of the total 
density matrix p(t) of the composite system is governed 
by some Hamiltonian of the form H — Hq + V, where Hq 
generates the free time evolution of the system and of the 
environment, and V describes the system-environment 
coupling. We work in the interaction representation and 
write the von Neumann equation of the combined system 
as 

j t p{t) = -i[V{t),p{t)]=C{t)p{t), (2.1) 

where V(t) is the Hamiltonian in the interaction picture, 
and C(t) denotes the corresponding Liouville superoper- 
ator. 

The projection operator techniques are based on the 
introduction of a projection superoperator V. This is a 
linear map 

p^Pp (2.2) 

which takes any operator p on the total state space H to 
an operator V p on H, and which has the property of a 
projection operator: 

V 2 = V. (2.3) 

Given a map V with this property one employs the pro- 
jection operator techniques to derive from the von Neu- 
mann equation (|2.1|) for the total density matrix p(t) 
exact and closed equations of motion for its projection 
Vp{t) (see Sec. IIIB|) . The basic idea underlying this ap- 
proach is the following. With an appropriate choice for 
the projection superoperator one intends to obtain a de- 
scription of the dynamics of system states which is much 
simpler and much more efficient than the description by 
means of the full density matrix p(t). Thus, the map V 
expresses the transition from the full representation in 
terms of the total density matrix p to a simplified, ef- 
fective description through a reduced set of dynamical 
variables defined by the structure of the projection V p. 

Equation (|2.3|) is our first basic condition. It is this 
condition which allows the formal application of the pro- 
jection operator techniques to open quantum systems. 
The ultimate goal is to determine form the equations 
of motion for V p(t) the dynamics of the density matrix 
ps(t) of the reduced open quantum system. To this end, 
we need one further condition. Namely, whatever the 
form of the projection superoperator is, we demand that 
V p contains the complete information required to recon- 
struct ps- We therefore impose the second basic condi- 
tion: 

p s = tr E p = tr E {Vp} . (2.4) 



The first equation is just the definition of the reduced 
density matrix which is obtained by taking the partial 
trace over the environment. The second equation states 
that, in order to determine ps, we do not really need the 
full density matrix of the total system, but only its pro- 
jection Vp. Hence, the reduced density matrix ps(t) is 
found by taking the environmental trace of the equations 
of motion for V pit). 

Within the standard approach to the dynamics of open 
systems using projection operator techniques one defines 
a projection superoperator of the form: 

Vp= {tr E p)®p E , (2.5) 

where pe is some fixed environmental density matrix 
referred to as the reference state. This superoperator 
clearly satisfies our basic conditions l|2.3[) and l|2.4|) . By 
use of this V the total state of the system is represented 
by means of the tensor product state ps ® Pe- Regard- 
ing the reference state pe as fixed, one uses the reduced 
density matrix ps(t) as the dynamical variable. Applying 
the projection operator technique one then finds a master 
equation for the reduced density matrix ps whose coef- 
ficients are given by certain multitime correlation func- 
tions defined by averages with respect to the reference 
state pe- In particular, the master equation obtained 
in second order of the coupling yields the Born approxi- 
mation of the dynamics which involves certain two-time 
environmental correlation functions. 

Our new class of correlated projection superoperators 
is obtained as follows. We take any orthogonal decom- 
position of the unit operator Ie on the state space of the 
environment, i. e., a collection of projection operators n a 
on He which satisfy 

n Q n fc = s ab n b , ^ n a = i E . (2.6) 

a 

Then we can define a linear map by means of 

Pp = ^tr B {n aP }®i-n a , (2.7) 

a 

where N a = tr^ {II a }. It is again easy to verify that this 
superoperator fulfills the requirements (|2.3|l and (|2.4|l . 
By contrast to the standard projection l|2.5|l which uses 
a representation by a tensor product state, the new pro- 
jection H2.7f) employs a set of un-normalized density ma- 
trices 

P f = ti E {n aP } (2.8) 

in order to describe the states of the composite system. 
The set of the density matrices Pg(t) therefore repre- 
sents the dynamical variables defined by the projection 
superoperator (|2.7|l . Applying the projection operator 
technique one is then led to a coupled system of equa- 
tions of motion for the Pg(t), from which one obtains 
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the reduced density matrix itself by means of the relation 
P s(t) = tr E {Vp(t)} = J2p ( s ) (t)- (2-9) 

a 

In the theory of entanglement (see, e. g., Ref. 01) 
a state of the form given b y E q. I|2.7I) is called separa- 
ble or classically correlated [43j . The approach based on 
a projection of this form thus tries to approximate the 
total system's state through a classically correlated but 
non-factorizing state. Of course, one can also construct 
projection superoperators which lead to non-separable 
(entangled) states Vp. The examples discussed below 
belong to the classes of projection operators defined by 
Eqs. (E3 and (|2~7|) . 

B. Equations of motion 

Basically, there are two variants of the projection oper- 
ator technique. The first one is the prominent Nakajima- 
Zwanzig method [THllTfil |. It leads to a first-order integro- 
differential equation for Vp{t) which contains a time 
integration over the past system history involving a 
certain memory kernel. The second variant is known 
as time-convolutionless (TCL) projection operator tech- 
nique [Tfll |23| which yields a time-local equation of mo- 
tion for Pp(t). We shall use this second variant of the 
projection operator technique in the present paper. Its 
advantage is that in any order of the coupling one only 
has to solve a first-order differential equation which is lo- 
cal in time. It should be emphasized, however, that the 
general considerations developed here may also be ap- 
plied to the Nakajima-Zwanzig projection operator tech- 
nique. 

The TCL projection operator method leads to an equa- 
tion of motion for the projection Vp(t) which is of the 
general form: 

j t Vp{t) = IC(t)P P (t)+I(t)Qp(0). (2.10) 

This is an exact inhomogeneous linear differential equa- 
tion of first order. Both the TCL generator IC(t) of 
the linear part and the inhomogeneity X(t) are explic- 
itly time-dependent superoperators. 

The inhomogeneous part of Eq. H2.10JI is determined 
by the projection Qp(0) of the initial state, where Q = 
I — V is the projection superoperator complementary to 
V, and / denotes the unit map. We observe that the 
inhomogeneous term vanishes if the initial state satisfies 
the relation Qp(0) — 0, i. e., if 

Vp(0)=p(Q). (2.11) 

Obviously, this relation simplifies the equation of motion 
(|2.1U|) . In the case of the standard projection (|2.5ll it 
implies that p(0) is an un-correlated tensor product state, 
i. e. p(0) = ps{0)®PE- In the case of the projection 12. 711 . 



however, Eq. (|2.11(l merely implies that the initial state 
is of the correlated form given by the structure of the 
projection superoperator (|2.7|) . 

In general, the TCL generator IC(t) and the inhomo- 
geneity T(t) are extremely complicated objects, and the 
exact solution of Eq. (|2.1U|I is as difficult as the solution 
of the full von Neumann equation for the total system. 
However, Eq. H2.10[l can be used as a starting point of 
a systematic perturbation expansion with respect to the 
strength of the interaction Hamiltonian V. With the help 
of the TCL technique one derives a closed expression for 
the corresponding expansion of the TCL generator: 

00 

K(t)=J2^n{t)- (2.12) 

n=l 

The n-th-order contribution is given by 

/■* rti rt n -2 

fcn(t) = dti dt 2 . ■■ dt n -i 

Jo JO Jo 

x (£(t)£(i!)£(t 2 ) . . . £(t„_i)) oc . (2.13) 

The quantities 

(£(t)£{t 1 )£(t 2 )...£{t n - 1 )) oc = 

^2(-l)«V£(t) . . . . . . £{t k )V ...V 

arc known as ordered cumulants [3(3, |32, 03 an d are 
defined by the following rules: (1) Write a string of the 
form V £ . . . CP with n factors of £ in between two Vs. 
(2) Insert an arbitrary number q of factors V between 
the £'s such that at least one £ stands between two suc- 
cessive Vs. The resulting expression is multiplied by a 
factor (— l) q and all £'s are furnished with a time ar- 
gument: The first time argument is always t. The re- 
maining £'s carry any permutation of the time arguments 
ti,t2,..., i„-i with the only restriction that the time ar- 
guments in between two successive 7-"s must be ordered 
chronologically. In the above expression we thus have 
t > • • ■ > ti, tj > . . . > tfc, etc. (3) Finally, the ordered 
cumulant is obtained by a summation over all possible 
insertions of V factors and over all allowed distributions 
of the time arguments. 

In many physical applications it may be assumed that 
the relations 

V£(t)£(h)...£{t 2n )V = (2.14) 

hold, which means that any string containing an odd 
number of £'s between successive factors of V vanishes. 
Following the above rules one then finds that all odd- 
order contributions lC2 n +i(t) vanish, while the second 
and the fourth-order terms take the form: 




dtxP£{t)C{h)V, (2.15) 
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and 

K 4 (t) = dtx dt 2 / dt 3 
Jo Jq Jo 

-vc(t)c(t 2 )vc(h)c(t 3 )r 

-VC{t)C{U)VC(t x )C{t 2 )V). (2.16) 

The performance of the formal expansion outlined 
above strongly depends, of course, on the choice of the 
projection superoperator V . In other words, the quality 
of the approximation obtained by truncating the expan- 
sion at a given order n crucially depends on the struc- 
ture of the chosen projection. It is important to note 
that the technique yields an expansion of a certain sys- 
tem of equations of motion, and not of the reduced sys- 
tem's density matrix itself. Taking different projection 
superoperators one uses entirely different sets of dynam- 
ical variables which obey completely different equations 
of motion. Hence, changing the projection superoperator 
amounts to changing the set of dynamical variables and 
the whole structure of the equations of motion, and to 
a non-perturbative re-organization of the expansion. It 
may even happen that the solution of the equations of 
motion in a given order for one particular projection rep- 
resents the solution to all orders for another projection. 
This point will be illustrated in Sec. IIVI bv means of a 
specific system-reservoir model. 

III. HILBERT SPACE AVERAGING 
APPROACH TO THE REDUCED DYNAMICS 

A. The Hilbert space average method 

The Hilbert space average method (HAM) is in essence 
a technique to produce guesses for the values of quan- 
tities defined as functions of a wavefunction \ip) if \ip) 
itself is not known in full detail, only some features of 
it. In particular it produces a guess for the expectation 
value (ip\A\ip) if the only information about is a set 
of different expectation values {tp\B n \tp) = B n . Such a 
statement naturally has to be a guess since there are in 
general many different \ip) that are in accord with the 
given set of B n but produce possibly different values for 
The question now is whether the distribution of 
{ip\A\ipys produced by the respective set of |^)'s is broad 
or whether almost all those \ip)'s yield {ip\A\ipYs that are 
approximately equal. It turns out that if the spectral 
width of A is not too large and A is high-dimensional 
almost all individual \ip) yield an expectation value close 
to the mean of the distribution of ( , 0|^4|'0)' s - In spite of 
this being crucial for the following we refer the reader to 
for details. To find that mean one has to average with 
respect to the \ip)'s. We call this a Hilbert space average 



A and denote it as 

A=[(^\Am m6nW=Bn y (3.1) 

This expression stands for the average of {tp\A\ip) over 
all that feature (ip\B n \ip) = B n but are uniformly 
distributed otherwise. Uniformly distributed means in- 
variant with respect to all unitary transformations that 
leave (ip\B n \%jj) = B n unchanged. One may rewrite l|3.1|) 
as 

A = tr{ia} with &=M{i>\} {W B nm=Bn y ( 3 - 2 ) 

How is a to be computed? Any unitary transformation 
that leaves {ip\B n \ip} — B n invariant has to leave a in- 
variant, i. e.: 

e td ae~ ld = a with [G, B n ] = 0. (3.3) 

This, however, can only be fulfilled if [G, a] = and this 
leads to the general form 

a = J2 b nB n . (3.4) 

n 

(In principal there could be addends of the form, e. g., 
B n B m , etc., but since all B n we are going to consider 
below together with zero form a group, those addends 
are already contained in the above sum.) 

Furthermore one of course has the following conditions: 

tv{aB m } = B m . (3.5) 

One thus obtains 

B m = J2tr{B m B n }b n (3.6) 

n 

from which the b n may be determined. Thus, the con- 
struction of a given Hilbert space average is defined with 
the help of Eqs. (|3.6p . lpT3)l . and According to this 

scheme "best guesses" for certain expectation values will 
be produced below. The explanation of HAM in full de- 
tail is beyond the scope of this text and can be found in 
00. 

B. HAM, projection operators and dynamics 

In the following we explain how HAM can be used 
to produce the reduced dynamics of a quantum system 
coupled to some environment, just like the techniques 
described in Sect. El Consider the full system's pure 
state at some time t, \^p(t)). Let D(t) be a time evolution 
operator describing the evolution of the system for a short 
time, i. e., \ip(t + r)) = D(r)|^(i)). This allows for the 
computation of a set of observables B n at time t + r: 

B n (t + t) = (iJj(t)\D\ T )B n D(T)\i;(t)). (3.7) 
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Now assume that rather than \t^(t)) itself only the set of 
expectation values B n (t) = (tfj(t)\B n \vp(t)) is known. The 
application of HAM produces a guess for the B n (t + r) 
based on the B n (t): 



B n (t+r) w {{4>\i)HT)B n D(i 



{(4>\B n \4>)=B n {t)y 



(3.8) 



(Note that here the \<j>) appear rather than the be- 
cause those are not actually realized states but denote 
the set of states over which the Hilbert space average 
has to be taken.) Iterating this scheme, i. e., taking the 
B n (t + t) for the B n {t) of the next step allows for the 
stepwise computation of the evolution of the B n 's. If the 
set of the B n is chosen such that it determines the local 
state of the considered quantum system completely this 
technique produces the local reduced dynamics. The re- 
sult is of course just like HAM itself only a best guess, 
but for appropriate systems this guess can be rather ac- 
curate. 

Here, the B n 's are chosen specifically as operators cor- 
responding to elements of the reduced density matrix of 
the considered system and the occupation probability of 
"energy bands" of the environment: 



B„ = Bj 



(3.9) 



where \i), \ j) are energy eigenstates of the considered sys- 
tem and n a is as described in Sect. Ill Al a projector, pro- 
jecting out the energy eigenstates of the environment be- 
longing to an interval AE a around some mean band en- 
ergy E a labelled by the index a. Let {tp\Bij a \i^) = Bij a , 
then one gets for the elements pij of the reduced density 
matrix: 



Pij 



- ^2 B 3ia- 



(3.10) 



Thus, the above defined set of expectation values deter- 
mines the reduced state of the system completely. The 
set of states a which belong to the Hilbert space average 
defined by the B^ a [in the sense of Eq. 1|3.2|) ] is, following 
the scheme described in Sect. IHlXl found to be 



N a ija ' 

yes 



(3.11) 



(This turns out to be the same state one would have 
gotten from minimizing the purity under the subsidiary 
condition set by the given expectation valu es Bj j a .) A 
comparison with the considerations of Sec. Ill Al reveals 
that a has exactly the form produced by the applica- 
tion of the projection superoperator V [see Eq. (|2.7(l ]. 
Exploiting Eqs. I|3.1|) . (|3.2|) . and 1)3. 11|) one can write a 
specific form of Eq. I|3.8[) for this case: 

B«a(* + T)» (3.12) 

]T ^tr{|m)(Z|n bJ Dt (r )| i)(j |n Q ^( r )} B lmb (t). 

Irnb 



Working in the interaction picture the dynamics of the 
full system is controlled by the interaction V(t). The 
time evolution is generated by the corresponding Dyson 
series. Thus, assuming weak interactions Eq. H3.12f) may 
be evaluated to second order in the interaction strength 
using an appropriately truncated Dyson series for D(t). 
This yields after extensive but rather straightforward cal- 
culations for the expectation values corresponding to di- 
agonal elements: 



B iia (t + T) = 

jb 



(3.13) 



Bjjbji) _ B iia (t) 
N b N a 



and for the expectation values corresponding to off- 
diagonal elements: 



B ija (t + T) 
Bij a (t) - 



(3.14) 



l B ija (t) 
2 N a 



^2(f(ikab,T)+f(kjab,r)), 



kb 



where 



f(ijab, t) 



2 dr' dr" 
Jo Jo 

xtv{Il a (i\V(T")\j)Il b (j\V(0)\i)} 



(3.15) 



Those f's are essentially integrals over the same envi- 
ronmental temporal correlation functions that appear in 
the memory kernels of standard projection operator tech- 
niques. But here they explicitly correspond to tran- 
sitions between different energy subspaces of the envi- 
ronment. (In Eqs. (|3.13() and l|3.14|l we assumed that 
correlation functions vanish unless they refer to correla- 
tions between parts of the interaction that are adjoints 
of each other as in Eq. I|3.15l) . We furthermore assumed 
tr{H a (i\V(T)\j)} — 0. Both conditions are not neces- 
sarily fulfilled but apply to the concrete model analyzed 
below.) Those correlation functions typically feature 
(short) decay times, i. e, integrating them twice yields 
functions which increase linear in time after the corre- 
sponding decay time. Thus, for r larger than the decay 
time Eq. (|3.15ll may be written as 



f(ijab,r) w N b ^(ijab)r, 



(3.16) 



where "/(ijab) has to be computed from Eq. I|3.15|l but 
typically corresponds to a transition rate as obtained 
from Fermi's Golden Rule. Especially it will only be non- 
zero for Ei — Ej « E a — E b for otherwise the correlation 
function oscillates rapidly before it decays and hence the 
corresponding integrals vanish. 

Inserting Eq. l|TTE|) into Eqs. (|3~T3|) and lITTljl and as- 
suming that the decay times of the correlation functions 
are small compared to the resulting decay times of the 
system (which are of the order of l/"f(ijab)), one can 
transform the iteration scheme into a set of differential 
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equations: 



model is taken to be H — Hq + V, where 



Ss 

H = A^cr x +Y]-^ni|ni)(m| 

— iVi 



dt 



B ija = - ^Bija ^ (l(kiba) + -f(jkba)) . 



(3.17) 
(3.18) 



A-fj 



This set of differential equations obviously determines the 
reduced dynamics of the considered system. It produces 
an exponential decay to an equilibrium state. Again 
those dynamics are only a guess, but as a guess they 
are valid for any initial state regardless of whether it 
is pure, correlated, entangled, etc. In contrast to the 
standard Nakajima-Zwanzig and TCL methods where 
initial states generally produce an inhomogencity [see 
Eq. l|2TTU)) ] which may be difficult to handle, HAM al- 
lows for a direct guess on the typical behavior of the 
system. However, a crucial condition for the application 
of the above scheme is that the decay times of the corre- 
lations are short enough such that even for larger times 
the evolution is well described by a Dyson series trun- 
cated at second order. This means that the scheme will 
break down altogether if the interaction is too strong. 



IV. APPLICATION 
A. The model 

To illustrate the general considerations of the previous 
sections we investigate the model of a two-state system S 
which is coupled to an environment E 9]. The environ- 
ment consists of a large number of energy levels arranged 
in two energy bands of the same width Se. The levels of 
each band are equidistant. The lower energy band con- 
tains Ni levels, the upper band iV 2 levels. The transition 
of the two-state system is in resonance with the distance 
AE between the bands (see Fig. . 



AE 



V 



system 



= 5e 

T 

i 

11 Se 

T 

environment 



iVi 



FIG. 1: a two-state system coupled to an environment con- 
sisting of two energy bands with a finite number of levels. 

The total Schrodinger picture Hamiltonian of the 



Se 



J2[&E+—n 2 )\n 2 ){n 2 \ 



(4.1) 



and 



V = \ c(n 1 ,n 2 )a + \n 1 )(n 2 \+h..c, (4.2) 

n 1 ,n 2 

Here and in the following the index n\ labels the levels 
of the lower energy band and n 2 the levels of the up- 
per band. a z and o+ are standard Pauli matrices. The 
overall strength of the interaction is parameterized by 
the constant A. The coupling constants c{n\^n%) are in- 
dependent and identically distributed complex Gaussian 
random variables satisfying: 



(c(ni,n 2 )) = 0, 
(c(n 1 ,n 2 )c(n' 1 ,n' 2 )) = 0, 
(c(ni,n 2 )c*(n' 1 ,ri2)) = S^^S^ 



(4.3) 
(4.4) 
(4.5) 



Transforming to the interaction picture we get the von 
Neumann equation H'2.1f) with the interaction picture 
Hamiltonian: 



V{t) = <r+B{t) + a-B*(t), 



(4.6) 



where 



B(t) =\J2 c(ni,n 2 )e-^^' n2 )*|n 1 )(n 2 | (4.7) 

"1,«2 



and 



w(ni, n 2 ) = Se 



n 2 ni 



N 2 



Ni 



(4.8) 



B. The standard approach 

1. Projection superoperator 

In the standard approach one uses a projection super- 
operator of the form given by Eq. I|2.5I) . Let us denote 
the projection onto the lower (upper) band by Hi (H 2 ): 



m 

H 2 = 5> 2 }<n 2 |. 



(4.9) 
(4.10) 



We consider initial states for which only the lower band 
is occupied: p(0) — ps(0)®Hi/Ni. Hence, if we take the 
reference state 



(4.11) 
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we have 



and 



Vp = (tr E p) ®p E =p s <g, —IIi, 

iVi 



Vp(0) = p(Q). 



(4.12) 



(4.13) 



In the following we write the elements of the reduced 
density matrix as 



Pij(t) = (i\ps(t)\j), 



hJ = 0, 1, 



(4.14) 



where |0) and |1) denote the lower and the upper level 
of the two-state system, respectively. If follows from 
Eq. I|4.13[) that the inhomogeneous term of the TCL mas- 
ter equation l|2.10[l vanishes. It can also be verified eas- 
ily with the help of the above forms for the projection 
superoperator and the interaction Hamiltonian that the 
condition i|2.14[l holds true. Thus, the TCL generator 
IC(t) contains only the contributions from even orders of 
the coupling strength A. 



2. TCL master equation of second order 

From the expression (|2.15|l for the second-order con- 
tribution of the TCL generator we find: 



K, 2 {t)Vp{t) = / <fti/ a (t,ti) 
Jo 

x [2v-Ps(t)v+ - W+<r-,Ps{t)}] ® Pe, 
where {•, •} denotes the anticommutator and 

/ 2 (Mi) = (tT E {B(t)B^(ti) PE }) 

= T*h(t-ti) (4.15) 

is the two-point environmental correlation function with 



72 



2ti\ 2 N 2 
Se 



(4.16) 



The angular brackets in Eq. I|4.15|l denote the average 
over the random couplings c{n\ 1 n2) which is determined 
by use of the relations H4.3fl - 14.5fl . The function h(r) 
introduced in Eq. I|4.15[) is then found to be 



h(r) 



Se sin 2 (fe-r/2) 
2^ (fc-r/2) 2 ' 



(4.17) 



where we have assumed a constant finite density of states 
for the environmental energy bands. This function ex- 
hibits a sharp peak of width 8s at r = and may 
be approximated by a delta function for times t which 
are large compared to the inverse band width, i. e., for 
fe ■ t > 1 we may approximate: 



h{tM) «72<5(*-ti). 



(4.18) 



This yields the second-order master equation for the re- 
duced density matrix: 



dt 



Ps(t) = 72 



<T-ps(t)<r+ - ^{cT + a-,ps(t)} 



(4.19) 



This is a quantum Markovian master equation in Lind- 
blad form, where the quantity 72 represents the Marko- 
vian relaxation rate. 

On the ground of the second order approximation one 
could naively expect that the master equation (|4.19(l pro- 
vides a reasonable approximation of the reduced system's 
dynamics if the relaxation rate 72 is small compared to 
the band width: 



72 < Se. 



(4.20) 



However, we are going to demonstrate that this is not 
true. A comparison with numerical simulations of the full 
Schrodinger equation and with the prediction of HAM 
shows that the long time dynamics is not correctly re- 
produced by this master equation. 



HAM 
Schrodinger 




3000 



FIG. 2: comparison of the numerical solution of the 
Schrodinger equation with the approximations given by HAM 
[Eq. (|4 . 2 1 p> ] and by the second and the fourth order of the 
standard TCL expansion [Eqs. (14.2311 and 14.311 1. Parame- 
ters: Ni=N 2 = 500, Se = 0.5, and A = 5 • 10" 4 . 

The approximation of HAM leads to the following ex- 
pression for the population of the upper level: 



Pii(*)=pn(0) 



7i 



71 +72 7i + 72 



' 2 c -(7l+72)« 



(4.21) 



where 72 is defined by Eq. H4.1bfl and we have introduced 
a further relaxation rate: 



7i 



2ttA 2 Ai 

sF~ 



(4.22) 



On the other hand, the TCL master equation H4.19(l gives: 

Pu(t)=pii(0)e-T 2t . (4.23) 

Thus, the TCL master equation predicts an exponential 
relaxation of the populations to zero, while the solution 
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obtained by means of HAM approaches the stationary 
population 



be: 



stat 
Pll 



Pn(0) 



7i 



7i +72 



(4.24) 



To judge the quality of the various approximations 
we have performed numerical solutions of the full 
Schrodinger equation corresponding to the Hamiltonian 
defined by Eqs. 14.111 and 14.211 . The initial state has been 
taken to be of the form |1) ® \x), where the environmental 
state \x) represents a superposition of the states \ni) of 
the lower band with independent Gaussian distributed 
random amplitudes of zero mean and equal variances. 
For certain parameter ranges we find an excellent agree- 
ment of the HAM prediction with the simulation results. 
An example is shown in Fig. [21 Note that for the param- 
eters of this figure we have 72 /Se = 3T0~ 3 , such that the 
standard Markov condition l|4.20|l is very well satisfied. 

We conclude that although the standard Markov con- 
dition (|4.20() is fulfilled the Markovian master equation 
14.191) does not yield a good approximation of the dy- 
namics for intermediate and long times. In particular, 
its prediction for the stationary state is totally wrong. 
The important point to note is that, in order to judge 
the quality of the Markov approximation, an analysis of 
the contributions of higher orders is indispensable. We 
also note that the same problem occurs if one uses the 
Nakajima-Zwanzig master equation. 



3. The master equation of fourth order and failure of the 
Born-Markov approximation 

To understand the failure of the Born-Markov approx- 
imation we investigate the fourth order of the TCL ex- 
pansion. The contribution of fourth order to the TCL 
generator is given by Eq. (|2.16|) . One finds that this 
contribution is determined by the two-point correlation 
function l|4.15|l and by the four-point correlation func- 
tion: 

hitMMM) = (tr E {S(t) J Bt(t 1 )S(i 2 ) J Bt(t 3 )p E }) . 

(4.25) 

The analysis shows that this function has two sharp peaks 
of width fe -1 at t = t\, t 2 = t$ and at t\ — £2, t = £3, 
and may be approximated, under the conditions of the 
previous section, by the expression: 

h{tMMM) « 72<*(* - *i)<*(*2 -*s) 

+7i 72 5(t - t 3 )<f(*i - i 2 ), (4.26) 

where 71 12 are defined by Eqs. (|4.22l) and (|4.1(jl) . 

The double-peak structure of the four-point correlation 
expressed by Eq. (|4.26|) has decisive consequences. With 
the help of the correlation functions given above the mas- 
ter equation of fourth order in the coupling is found to 



d , . 



T(t) 
f(t) 



cr-ps(7+ - tjW '-' As} 



(4.27) 



(T+<T-ps<7+<J- - -{<7 + cr_, ps} 



where 



r(t) = 72(1-71*), f(i) = 7i72t 



(4.28) 



Equation l|4.27|l is a master equation with time- 
dependent relaxation rates T(t) and T(t). The second- 
order contribution to the rate T(t) is given by 72, while 
the fourth order yields the contribution 72 • (71 1). There- 
fore, the fourth-order term of the TCL generator is small 
compared to the second-order term only if the times t 
considered satisfy the additional condition 



Tit < 1. 



(4.29) 



Thus we see that the occurrence of terms proportional 
to t is responsible for strong deviations from the Marko- 
vian behavior. These terms are due to the double-peak 
structure of the four-point correlation function f&. We 
remark that this structure is markedly different to the 
usual situation of the coupling of an open system to a 
Bosonic field vacuum, for example. In this case fn has 
only a single peak and, hence, the above phenomenon of 
strong deviations from the Markovian dynamics for weak 
couplings does not occur. 

The master equation l|4.27|l yields the coherences: 



Poi(t) =poi(0)e^ 2t/2 , 
and the populations: 

pii(t)=p n (0)e-^ t+ ^ t2 / 2 . 
For 71. 2 1 -C 1 we find the expansion: 



(4.30) 



(4.31) 



Pii(t) =Pn(0) 



1 9 
1 - 72* + ^72(71 + 72 )t + 



which is seen to coincide with the corresponding short- 
time expansion of the HAM approximation given by 
Eq. H4.21f) . Thus we conclude that the TCL expansion 
based on the standard projection reproduces the short- 
time behavior predicted by HAM within the given orders. 
Correspondingly, the TCL approximation of fourth order 
clearly improves the approximation for short times, but 
leads to un-physical results for longer times and diverges 
in the limit t — > 00, as is illustrated in Fig. [21 We note 
that a similar situation occurs for the spin star model 
studied in Ref. |45| . 

Summarizing, the fourth order clearly indicates that 
the TCL expansion does not converge uniformly in t. It 
only provides a short-time expansion of the dynamics. 
As a result of the emergence of terms which are given by 
powers of t, it is impossible to obtain valid predictions on 
the long-time dynamics if one truncates the TCL series 
at any finite order. 
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C. TCL expansion using the correlated projection 
superoperator 

In view of the above analysis the following question 
arises: Is it possible to construct a new projection super- 
operator V whose corresponding TCL expansion yields 
the full prediction of HAM already in lowest order, and 
leads to a systematic expansion around HAM in higher 
orders? To answer this question we consider the following 
projection superoperator: 



Vp 



tv E {nip} «> -^ni + tr B {u 2 p} <s> ^-n 2 

iVl 1\ 2 



(4.32) 



This projection belongs to the class of superoperators in- 
troduced in Eq. (|2.7|) . By this ansatz the total system's 
state is approximated by a separable but non-factorizing 
state. The dynamical variables are the un-normalized 
density matrices p^ and p^ which are correlated with 
the projections onto the lower and the upper band, re- 
spectively. The reduced density matrix of the two-state 
system is found with the help of Eq. I|2.9I) : 



Ps (*) = tv E { V P (t)} = pW (t) + (t) 



,( 2 ) 



(4.33) 



Assuming that the initial state is of the correlated form 

p(0) - p { s ] (0) ® ^-Hx + p<g\o) ® in 2 , (4.34) 

we have VpiQi) — p(0) and the inhomogeneous term of 
the TCL equation l|2.10|l vanishes. 



1. The second-order master equation 

Using the new projection superoperator (|4.32|) we get 
the following second-order TCL equation: 

j t v P (t) = p %\t)®±n 1 + pP(t)®±-ii 2 

= JC 2 {t)V P {t), (4.35) 
where the TCL generator takes the form: 

JC 2 {t)V P {t) = (4.36) 



dtxh(t - h) 



+ / dtxhit-h) 



2~i 2 a-p [ g ] a + - 7i{er_cr+,p^} 



(2)l 



iV 2 



-n 2 



Combining Eqs. (|4.35(l and (|4.36|l and assuming again 
fe • f > 1 we find the equations of motion: 



dt 



d p ( s 2) (t) 



dt 



(2) 

(1) 

l 2 V-Ps'<J + 



21 
2 

7i 
2 



{<r +C r_,p«}, (4.37) 
{a_a + ,pf}. (4.38) 



This is a coupled system of first-order differential equa- 
tions for the two density matrices p^g\i) and Pg(t). The 
elements of these matrices are written as 



p^(t) = m\t)\j)> 



p^v 



- /.•!„« 



i\pf{t)\j). (4.39) 



Schrodinger 
new TCL2 
HAM 




100 200 300 400 

t 



FIG. 3: comparison of the second-order TCL approximation 
using the new projection superoperator [Eq. 14.511 1. of the 
approximation given by HAM [Eq. 14.21H 1. and of the nu- 
merical solution of the Schrodinger equation. Parameters: 
^1=^2 = 500, 6e = 0.5, and A = 0.001. 

The equations (j4.37|l and lj4.38J) can now be used to de- 
rive an equation of motion for the reduced density matrix, 
making use of Eq. (|4.33|l . First, we get from Eq. I4.38|l : 



d 



(4.40) 



Jt P$(t) = 12^(^-^(1). (4.41) 

We assume again that initially only the lower band is 
populated: 



,4 2) (o) = o. 

It thus follows from Eq. (|4"4TJ|) that 



(4.42) 



(4.43) 



From Eq. I|4.37|l we find 



dt 



p[\\t)=JiP$(t)- 72 p[\\t). (4.44) 



From Eqs. 14.4411 and (|4.41|) we see that the quantity 

Putt) + Pro (*) i s constant. With the help of the initial 
condition (|4.42l) we thus have 

P$(t)=P$(0)-P$(t)- (4-45) 
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Substituting this into Eq. I|4.44[) we find: 



PS® 



-(7i+72>S ) (*)+7iPu(0)- (4-46) 



since pn (t) = p[\\t) because of Eq. i|P3|l . we finally 
arrive at the equation of motion for the populations: 



dt 



Pn(t) = -(71 + 72)pn(<) + 7iPii(0)- 



(4.47) 



In a similar manner one is led to the equation of motion 
for the coherences: 



d M ^2 M 
Jt Poi(t) = - Y Poi(t). 



(4.48) 



The dynamics of the reduced density matrix is thus de- 
termined by Eqs. Ij4.47jl and Ij4.48jl . These are time- 
local first-order differential equations with constant co- 
efficients. They are identical to the equations of mo- 
tion obtained using HAM. In particular, the solution of 
Eq. (|4.47|l is given by the expression l|4.21l) . Hence, we 
conclude that the lowest order of the TCL expansion with 
the projection superoperator introduced in Eq. I|4.32ll in- 
deed reproduces the HAM prediction. 




Schrodinger 
new TCL2 



HAM — - 



10 20 30 40 50 60 70 
t 

FIG. 4: the same as Fig.EJfor A = 0.003. 



Schrodinger 
new TCL2 
HAM 




FIG. 5: the same as Fig.EJfor A = 0.01. 

We observe that the dynamics of the populations p\\ (t) 
is strongly non-Markovian because of the presence of 



the initial condition pn(Q) on the right-hand side of 
Eq. (|4.47|) . This term expresses a pronounced memory 
effect, namely it implies that the dynamics of the popu- 
lations never forgets its initial data. Note also that the 
dynamics of the reduced density matrix is not in Lind- 
blad form and does even not represent a semigroup. It 
does, however, lead to a positive dynamical map, as can 
easily be verified. 

In the transition from Eq. I|4.36|l to the Eqs. (|4.37|) 
and (|4.38|) we have assumed for simplicity that the times 
t considered satisfy the condition 5e-t ^ 1. Without this 
condition the master equations H4.37[) and (|4.38|) must be 
replaced by 



/ dtih(t-ti) (4.49) 
Jo 



dhh(t - h) 



2 l2 a-pf{t)a + - 7l {<7_<7+,p^(i)} 



(4.50) 



(2), 



These equations describe the full time dependence as it is 
predicted by the TCL expansion in second order. They 
lead to the following populations of the upper level: 



Pii(*)=Pu(0) 



71 



72 



-r(t) 



71 +72 7i + 72 



(4.51) 



where 



T{t) =2(71+72) / dt x / dt 2 h(ti-t 2 ), (4.52) 
Jo Jo 

and the function h(r) is given by Eq. (|4.17|) . 

In Figs. |31 El and we compare the result given by 
Eq. 1)4. 51|) with the prediction of HAM and with numer- 
ical simulations of the Schrodinger equation. The figures 
clearly show that already the lowest order of the TCL 
expansion with the new projection superoperator gives a 
good approximation of the dynamics. It not only yields 
the correct stationary state, but also reasonable predic- 
tions on the relaxation times even for rather strong cou- 
plings, where deviations from the HAM result are large. 
Note that the rates 71 2 of Fig. |3 differ from those of 
Fig. [SI by a factor of 100, and that the parameters of 
Fig. El correspond to the ratio 7^2/fe ~ 1.3. 



2. The master equation of fourth order 

We have seen that the lowest order of the TCL expan- 
sion obtained with the help of the correlated projection 
superoperator is capable of reproducing the prediction 
of the HAM approximation, and even improves this ap- 
proximation for larger couplings. The question is now: 
What happens in higher orders of the expansion? As our 
numerical simulations indicate, higher-order corrections 
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should be small for all times. We show that, by contrast 
to the case of the standard projection, this is indeed the 
case. 



J 



The fourth-order contribution to the TCL generator is 
obtained by using the projection superoperator defined 
by Eq. 14.32(1 in the general expression ((2.16(1 . One finds: 



K+(t)Vp = 



t i>tl pt2 

dt\ I dt 2 / dt 3 
(i Jo Jo 



2(7? + 7172) [h(t - t 2 )h(h - t 3 ) + h(t - t 3 )h(h - t 2 )\ a+pfa- 
-2 7 i 72 [h(t - t 2 )h(ti - i a ) + 2h(t - t 3 )h(h - i 2 )] a+a-p^a+a- 
- [llKt - t 2 )h(tt - i 3 ) + ( 7 | - 7i72)/i(i - h)h{t x - t 2 )} {a + <j_,pf} 



2(72 + 7i72) [h(t - t 2 )h{h - t 3 ) + h{t - t 3 )h{h - t 2 )] tr-p^a 



-2 7l72 [h(t - t 2 )h(h - t 3 ) + 2h(t - t 3 )h{h - t 2 )] a-(j +P fa-(j + 
- [llKt - t 2 )h(h - i a ) + ( 7l 2 - lll2 )h{t - t 3 )h(h - i 2 )] {a-a+, p^} 

I 



-n 2 



Here, /i(t) is again a function which is sharply peaked 
at r = and may be replaced by the delta function 
<5(t) for times which are large compared to the inverse 
band width. The decisive point is the following. When 
carrying out the three-fold time integrations in the above 
expression, no terms emerge which grow like powers of t. 
This is due to the fact that the integrands do not contain 
terms of the form h(t 2 — t 3 ). As a result, the contribution 
from K.i{t) remains small for all times i, and the limit 
t — > 00 of the generator exists. 

To illustrate this point let us model the function /i(t) 

by 



where we have introduced the rates: 



Mr) = £e-*-H 



(4.53) 



which approaches S(t) for infinite band width. The time- 
integrations may then easily be carried out to give: 



J dh j 1 dt 2 JJ dt 3 h(t - tjhfa - i 3 ) pa 



dh / dt 2 / dt 3 h(t - t 3 )h(t! - t 2 ) 
Jo Jo 



1 

8fe' 



Hence we see that K.i(t) becomes time- independent for 
fe-i 3> 1 and that the fourth order of the expansion leads 
to the equations of motion: 



d Mm -n (2) r 2 r (1)-, 
jlPs (*) = IV+Ps <*- - y{0+0"-,P]s } 

(i) 



-r 3 a + a-p ( s 'a+(7-, 



(4.54) 



d (2)/ + \ p (!) r i ( (2)1 



-rw- 



(2), 



3C-cr + p s (T_(T +) 



(4.55) 



r 2 
r 3 
fi 

f 2 



71 
72 

f 3 
7i 

72 



7i +72 



2Se 

272 — 7i 



4& 
37i72 
4& ' 
271 - 72 



A5e 
71 + 72 



25e 



(4.56) 
(4.57) 
(4.58) 
(4.59) 

(4.60) 



This shows that the fourth order merely yields correc- 
tions of order 0(71,2/ 5s) to the equations of motion 
found in second order. Thus, as expected the TCL se- 
ries obtained with the correlated projection superopera- 
tor indeed leads to a systematic perturbation expansion 
around the approximation suggested by HAM. 

We finally mention that the equations of motion for the 
populations of the reduced system in fourth order take 
the form: 



dt 



Pll (t) = -(ri + f 2 )p n (i) + ripn(o), 



which is easily solved to yield: 



Pii(*) = pn(0) 



7i 



7i +72 7i + 72 



_22_ e -(r 1+ f 2 )t 



(4.61) 



(4.62) 



The interesting point to note is that the stationary state 
is identical to the one found in second order. Thus, the 
stationary state is not affected by the fourth-order cor- 
rections. However, the rate of the relaxation into the 
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stationary state is found to be 



I\ +f 2 = (71 +72) 



7i + 72 



2Se 



(4.63) 



which is seen to be larger than the relaxation obtained 
in second order. 



V. CONCLUSIONS 

In this paper we have analyzed non-Markovian quan- 
tum processes by means of the time-convolutionless pro- 
jection operator technique and of the Hilbert space av- 
erage method. It has been demonstrated that by use 
of a class of projection operators which project onto 
correlated system-environment states, an efficient non- 
perturbative treatment of certain non-Markovian pro- 
cesses is possible. The correlated projections have been 
shown to correspond to the idea of a best guess underly- 
ing the approximation of HAM. 

The general mathematical formalism of the projec- 
tion operator techniques does not tell us which projec- 
tion superoperator should be used for a given system- 
environment model. The choice of an appropriate pro- 
jection V depends on the structure of the model under 
study and must be based on physical considerations. The 
aim is of course an efficient description, i. e., a descrip- 
tion which can be expected to yield accurate results even 
in low orders of the coupling. Once a certain projection 
V has been chosen, one can use the perturbation expan- 
sion in order to check explicitly whether or not higher 
orders remain small and, thus, whether or not V enables 
a computationally efficient treatment of the reduced dy- 
namics. 

We emphasize that the projection operator techniques 
yield an expansion of the equations of motion for the 
relevant variables V 'pit), and not an expansion of these 
variables itself. Different projection superoperators lead 
to different sets of relevant variables and, hence, to 
equations of motion with completely different structures. 



Consequently, the usage of different projections implies a 
complete re-organization of the perturbation expansion. 
This point has been demonstrated here by means of a 
specific system-environment model. As we have seen, the 
TCL technique which is based on the correlated projec- 
tion superoperator yields accurate results for this model 
already in lowest order, while the standard TCL proce- 
dure fails in any finite order of the coupling. 

Our results suggest applications to other models show- 
ing strong non-Markovian effects. For example, it is clear 
that the class of correlated projections introduced here 
may also be applied to a generalization of the model 
studied in Sec. IIVI which involves any number of well- 
separated environmental energy bands. 

In Sec. Ill Al we have formulated two general conditions 
for suitable superoperators V. The first one [Eq. (|2.3[) ] 
is simply the condition that the map V be a projection 
operator. The second condition [Eq. (|2.4|) ] requires that 
the projection Vp of any state p of the composite sys- 
tem contains the full information which is necessary to 
extract the reduced density matrix p$ from the relevant 
variables. It is easy to construct classes of superopera- 
tors which satisfy these two conditions and are even more 
general than the class of correlated projections investi- 
gated here. As mentioned already the latter project the 
total state onto a separable, classically correlated system- 
environment state. This suggests exploiting the possibil- 
ity of using superoperators which project onto nonsepa- 
rable, entangled quantum states. In this way one might 
be able to investigate the dynamical significance of en- 
tanglement in non-Markovian quantum processes. 
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